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Abstract 


The 5 -ray theory of (rack structure is compared with experimental data for the 
radial dose from heavy ion irradiation. The effects of electron transmission and the 
angular dependence of secondary electron ejection are included in the calculations. 
Several empirical formulas for electron range and energy are compared in a wide 
variety of materials in order to extend the application of the track-structure theory. 
The model of Rudd for the secondary electron spectrum in proton collisions , which is 
based on a modified classical kinematics binary encounter model at high energies 
and a molecular promotion model at low energies , is employed. For heavier projec- 
tiles , the secondary electron spectrum is found by scaling the effective charge. Radial 
dose calculations for carbon , water , silicon , and gold are discussed. The theoretical 
data agreed well with the experimental data. 


Introduction 

The 5-ray theory of track structure attributes the 
radiation damage from and the detection of the passage 
of heavy ions through matter to the ejection of electrons 
(5 rays) from the material by the passing ion (refs. 1-5). 
The track-structure theory has a long history of providing 
the correct description of a variety of phenomena associ- 
ated with heavy ion irradiation. Track-structure theory 
provided the first description of the spatial distribution of 
energy deposition from ions by a formula for the radial 
distribution of dose, as introduced by Butts and Katz 
(ref. 1) and Kobetich and Katz (ref. 2). This description 
led to many experimental measurements of the radial 
dose (refs. 6-10). The response of physical detectors, 
such as organic scintillators (ref. 5), thermoluminescence 
detectors (TLD’s) (ref. 5), alanine (ref. 11), nuclear 
emulsion (ref. 12), and the Fricke dosimeter (refs. 5 
and 13), to heavy ions has been described with track- 
structure theory. Many biological effects, such as the 
thindown of mammalian cells (ref. 5), which was pre- 
dicted nearly 20 years prior to the first experimental 
observation (ref. 14) have also been described. Track- 
structure theory is used to develop improved lithography 
methods by using ion beams for applications in micro- 
electronics and microtechnology (ref. 15). 

The radial dose distribution and the geometry of a 
target site are used in track-structure theory to map y-ray 
response to ion response. The radial dose for intermedi- 
ate distances from the track structure is known to 
decrease with the inverse square of the radial distance to 
the path of the ion, which has led to simplified formulas 
for many applications (refs. 1,5, 16, and 17). The radial 
dose both near and far from the path of the ion is difficult 
to predict because of uncertainties in the electron range 
and energy relation, the angular dependence of the 
secondary electron production cross section, and the 
effects of 5-ray transport in matter, especially for 


condensed-phase matter. However, many track-structure 
calculations have used simple, analytic forms for the 
radial dose from ions. The electron transmission and the 
angular dependence of electron ejection are ignored and 
simplified electron range-energy relations are used. In 
this paper, these factors are considered by following the 
method described in references 2, 3, and 17 and compari- 
sons are made to experimental data for radial dose distri- 
butions. Substantial changes in the physical inputs of the 
calculations were made. These changes include an 
improved model for the secondary electron spectrum for 
proton collisions with atoms and molecules (ref. 19) and 
improved formulas for the electron range-energy and the 
stopping power. 

Radial Dose Formalism 

The calculation of the radial dose D (t) as a function 
of the radial distance of the path of an ion of charge num- 
ber Z and velocity P was formulated in references 2 to 4 
and 18. In formulating the spatial distribution of energy 
deposition as charged particles pass through matter, the 
dominant mode of energy loss is assumed to be ioniza- 
tion due to electron ejection from the atoms of the target 
material. Electrons of range r that penetrate into a mate- 
rial a distance / have residual energy W, which is given 
by the energy co to go the residual range r - 1 . The resid- 
ual energy of an ejected electron (5 ray) is written in 
functional form as 

W(r,t) = CD (r - t) (1) 

In equation (1), r is the practical range (determined by 
extrapolating the linear portion of the absorption curve to 
the abscissa) of an ejected electron with energy ca. When 
the range-energy relation in a given target material is 
known, the residual energy is then evaluated with 
equation (1). 



The energy dissipated E at a depth t by a beam con- 
taining one electron per cm 2 is represented in reference 2 
as 


E = | (11 W) (2) 

where T] is the probability of transmission of the 
electrons. 

As noted in reference 4, equation (2) neglects sev- 
eral effects. First, it may neglect backscatter, although it 
may be argued that the energy lost from a layer dt by 
backscatter is compensated by energy gained from back- 
scatter from later layers. Second, all' electrons are repre- 
sented by an underscatter class. Third, the energy 
deposited by the electrons that penetrate to a thickness 
t > r is neglected. Such shortcomings could be overcome 
by direct solution of the electron transport (ref. 20) or 
through the use of Monte Carlo methods (ref. 21). How- 
ever, the model of Kobetich and Katz from reference 2 
has the advantage of simplicity with reasonable 
accuracy. 

The transmission function used is based on the 
expressions of Dupouy et al. (ref. 22) as modified by 
Kobetich and Katz (ref. 4) and is given by 

Tl(r,0 = exp[-(< 7 //r) p ] (3) 

with 

q = 0.0059 Z°' 98 + 1.1 (4) 

and 

p = 1.8 (log !0 Z r )“ 1 + 0.31 (5) 


In equation (6), Z is the effective charge number of the 
ion, which was represented in reference 24 as 


Z 


* 


= z 


1 - exp 



( 8 ) 


The electron-binding effects were considered by 
Kobetich and Katz (ref. 2) after the experimental find- 
ings of Rudd et al. (ref. 25), who found that a> may be 
interpreted as the total energy imparted to the ejected 
electron with a kinetic energy of W , so that co in 
equation (6) is replaced by 

co = W + I (9) 

Results from equation (6) must be summed for compos- 
ite materials. The average charge and mass number and 
the density of several materials are listed in table 1. Val- 
ues of mean excitation 7/ from references 26 and 27 and 
values for electron density Nj are listed in table 2. 

Rudd has provided a parameterization of the electron 
spectrum after proton impact that was based on a binary 7 
encounter model modified to agree with the Bethe theory 
at high energies and with the molecular promotion model 
at low energies (ref. 19). For water, the contributions 
from five shells are included (ref. 19). This model of 
Rudd is considered herein and scaled to heavy ions by 
using effective charge. In figure 1, the secondary 7 elec- 
tron spectrum from equation (6) and the model of Rudd 
for several proton energies are shown. Large differences 
between the models occur for all electron energies below' 
proton energies of about 1 MeV and for small electron 
energies at all proton energies. 

Electrons of energy co are ejected at an angle 0 rela- 
tive to the path of a moving ion and described by classi- 
cal kinematics as 


where Z T is the atomic number of the target material, and 

1 7 

r and t are in units of g/cm . 

In order to estimate the number of free electrons 
ejected by an ion per unit length of ion path with energies 
between co and co + <7co , the formula given by Bradt and 
Peters (ref. 23) was used by Kobetich and Katz (ref. 2) 


dn _ 2nNZ* 2 g 4 1 

//co 2^2 2 
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(6) 


where e and m are the electron charge and mass, N is the 
number of free electrons per cm 2 in the target, and co /?j is 
the classical kinematic value for the maximum energy 
that an ion can transfer to a free electron and is given by 


co 


m 


2 mcV 
1-P 2 


(7) 


cos 2 e = — do) 

CO 

m 

for the collision between a free electron and an ion. 
Equation (10) indicates that close to the path of the ion, 
where distances are substantially less than the range of 
5 rays (co«co m ) , it is sufficient to consider that all 
5 rays are normally ejected. The energy that the electrons 
dissipate in cylindrical shells with an axis along the path 
of the ion may then be found from the energy dissipation 
of normally incident electrons. If the 5 rays far from the 
path of an ion have an important role in a particular 
response, then the angular dependence, as well as the 
dependence of electron range, on the velocity of the ion 
becomes crucial. 

If £ is the energy flux carried by 5 rays through a 
cylindrical surface of radius t whose axis is the path of 
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the ion, the energy density E deposited in a cylindrical about the classical kinematic ejection value described by 

shell of unit length and mean radius t is given by equation (10), so that 


-1 dt 
Intdt 


(ID 


5 ? 3 H -£«*“> 


(15) 


The total energy flux is found by integrating the energy 
flux carried by a single electron, given by r\W , over the 
distribution of the 5 ray and summing the contributions 
of all atoms in the material 

e(0 = Xj“ w ‘ / '^W ( r,co ) ti ( t,co)^ (12) 

i 1 

In equation (12), the lower limit co ; is the energy for an 
electron to travel a distance t, and the upper limit co^ - L 
is the maximum kinetic energy that can be given to the 
electron by the passing ion. Using equation (11) and 
equation (12), the energy density distribution may be 
written as 

i 

and Eft) is identified as the radial distribution of dose. 

To consider the angular dependence of the ejected 
electrons, the energy deposited by a ray ejected at an 
angle 0 in a cylindrical shell of radius t centered on the 
path of the ion is assumed to be the same as the energy 
deposited by an electron normally incident on a slab at 
depth f/sin0 , as shown in figure 2. Kobetich and Katz 
assume that differences between the geometry of the slab 
and the cylindrical shell do not greatly affect the energy 
density distribution, because the differences in the 
energy density at t caused by the electrons scattered in 
path A are compensated by those scattered in path B of 
figure 2 (ref. 3). 

The energy density distribution, which includes an 
angular distribution of the ejected electrons, is assumed 
to be 


where 

/(9,(0) = (16) 

[0-e (co) ] 

CO 

with Q c (to) determined as the root of equation (10), A is 
a normalization constant, and A" is a constant. The con- 
stant K may have some dependence on the energy of the 
incident ion and target material; however, K is estimated 
as 0.015 keV from the data of references 28 to 31. 
Results of equations (15) and (16), which use the model 
of Rudd (ref. 19) for dnidtti in equation (15), are shown 
in figures 3(a) and 3(b). 

Range and Energy Formula in Arbitrary 
Media 

The electron range and energy relationship is diffi- 
cult to evaluate theoretically and, because of the com- 
plexity of the electron transport problem, empirical 
expressions based on experimental measurements have 
been developed (refs. 2, 4, and 32-36). Over a limited 
energy range, a power law of the form r - £co a will be 
approximately correct and is used in references 1, 34, 
and 36. The residual range of the power law is easily 
found by inversion and leads to an analytic form for the 
radial distribution of dose with the simplifying assump- 
tions of normal ejection and unit electron transmission. A 
more accurate form, which is given in reference 32 and 
modified in reference 4, is the formula 

r = < I7) 


E(t) = 


-1 
2j it 


ZJ* J 


“""Vo | 

(0) dt 


x [T| (r, (0, 0) W (t, co, 0)] 


dn. 


d(0 dQ 


(14) 


where 

A = [o.81Z“°' 38 + 0.18]xlO“ 3 g(cm 2 kev) ' (18) 


The angular dependence of (0 {i 7], and W is shown in 
equation (10). Experimental measurements for the dou- 
ble differential cross section of electron ejection are 
available for only a few ions and mostly at modest ion 
energies of <10 MeV/amu (refs. 28-31). 

A qualitative model for the angular distribution of 
the secondary electrons assumes that distribution peaked 


B = 0.21 Z 7 -°' 555 + 0.78 (19) 

C - (l.iz5: 29 + 0.21 jx 10' 3 keV _1 (20) 

which was determined by extensive comparison with 
experimental data for practical range in many materials. 
Equation (17) is inverted to provide co = co (r) . 
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As a final parameterization, the range formula of 
reference 33 is considered 


r = a. 


i a 2 


log ( 1 + a 2 t) - a 3 x^ 1 + a 4 1 


where 


X = co/m 

b 2 

o. j — b^Aj/Zj- 


( 21 ) 


(ref. 37), Kobetich and Katz (ref. 4), and Tabata et al. 
(ref. 33) for water are also shown. In figure 8, dWidr 
from equation (17) or equation (21) is compared with 
experimental data from references 38 and 39 for stopping 
power of an electron in water. The model of Tabata et al. 
(ref. 33) agrees well with experimental data to about 
1.0 keV. This model will be used for the electron range 
energy in radial-dose calculations. 

Calculations of Radial Dose 


a 2 ~ ^3 
°4 = b 6~ b ! Z T 


a 5 ~ b%/Z T 


(22) 


The values of b, from reference 33 are listed in table 3. 
Equation (21) reduces to equation (17) when a 1 x« 1 
and a 5 = 1 . A parameterization of the inversion of 
equation (21) is provided in reference 33 as 


In figures 9(a) and 9(b), the effects of electron trans- 
mission on calculations of radial dose in water are 
shown. Calculations are for proton projectiles; however, 
the radial dose is determined approximately by Z 2 /| 3 2 
from which results for other ions can be found. 
Figures 9(a) and 9(b) illustrate that the transmission fac- 
tor affects the radial dose calculation only very close to 
and very far from the path of the ion. The normalization 
and the expected decrease of radial dose with increasing 
distance as 1 / 1 2 are unchanged by including the trans- 
mission factor. 


x = c, 


exp 


c 2 + cy^l+c 4 r 


(23) 


where 


Cj = d^/Zj 

dy/Aj 

c 2 - d 2 Z T 
c 3 = d 4 -d 5 Z T 


c 4 ~ d^/Z T 


c 5 ~ d^/Z T 


(24) 


with d { listed in table 3. 

A logarithm and polynomial relationship has been 
used by Iskef, et al. (ref. 35) and more recently by Zhang, 
et al. (ref. 36). This, however, is less useful for the radial 
dose model because the inversion formula for 
co = co (r) is not found easily. 

In figures 4 to 6, electron range and the energy of the 
ion from equation (7) is plotted for several materials. In 
the low to intermediate energy range, the formula agrees 
closely; however, large differences occur below 1 MeV/ 
amu, especially for lighter materials. Above 1000 MeV/ 
amu, large differences also occur which grow with the 
increasing charge of the target material. In figure 7, the 
formulas of Iskef et al. (ref. 35), Waligorski et al. 


In figure 10, the radial dose calculations are com- 
pared with experimental data from references 6 to 9 for 
several projectiles with ion energies from 0.25 to 
377 Me V/amu. Figure 10 illustrates the decrease in 
radial dose with increasing distance 1 /t 2 in the interme- 
diate distance range. Close to the ion track (t < 10 nm) 
a contribution to the radial dose from molecular excita- 
tions, as discussed in reference 37, is expected but is not 
included in the present calculations. It is important to 
keep the contributions from excitation and ionizations 
distinct, because the secondary electron dose from ion- 
ization is assumed to be responsible for physical effects 
by heavy ions. 

At large distances, the inclusion of angular depen- 
dence in equation (14) offers a substantial improvement 
in the accuracy of calculations. Equation (15) provides 
an improvement in accuracy over the ejection angle 
model of classical kinematics (eq. (10)) at lower energies 
(<2MeV/amu). At higher energies, equation (15) 
appears to underestimate the radial dose at large dis- 
tances. Clearly, more information on the double differen- 
tial cross section for electron ejection is required. 

In figures 11(a) to 11(c), the effects of the radial 
dose calculations for several velocities in carbon, silicon, 
and gold are illustrated. The data shown in figure 1 1 
were determined with equation (6) from the secondary 
electron spectrum. This model is capable of providing 
the radial dose for an arbitrary ion in a wide variety of 
materials. 
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Concluding Remarks 

A model for the radial distribution of energy depos- 
ited about the path of a heavy ion developed prior to 
most experimental measurements of this distribution was 
improved. Theoretical results from the improved model 
were compared with experimental data for a variety of 
ions. Improved models of electron-range energy and 
stopping power and the electron-ejection spectra and 
angular distribution were used in calculations. Excellent 
agreement with experimental data was found; however, 
more information on the double differential cross section 
for electron ejection is required. Calculations of the 
radial dose from heavy ions in several materials of inter- 
est for spacecraft design and microelectronics are pre- 
sented. The radial dose model developed in this report is 
useful in determining the response of many detectors and 
components to space radiations. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
December 2, 1994 
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Table L Average Atomic Number (Z), Average Mass Number (A), and Density of Materials 


Material 

Composition 3 

Z 

A 

Density, g/cm 3 

Carbon 

c 

6.0 

12.0 

1.95 

Aluminum 

A1 

13.0 

27.0 

2.702 

Silicon 

Si 

14.0 

28.0 

2.33 

Copper 

Cu 

29.0 

63.0 

8.92 

Tin 

Sn 

50.0 

120.0 

7.28 

Gold 

Au 

79.0 

197.0 

19.30 

Lead 

Pb 

82.0 

208.0 

11.34 

Lexan Polycarbonate 

C16H14O3 

6.1 

12.1 

1.20 

Cellulose Nitrate 

c 6 h 8 o 9 n 2 

7.02 

14.0 

1.35 

Water 

h 2 o 

7.22 

14.3 

1.0 

Air 

0.755N, 0.2320, 0.013Ar 

7.52 

15.1 

0.001293 

Quartz 

Si0 2 

10.8 

21.7 

2.66 

Muscovite Mica 

KAl 3 Si 3 0 16 (OH) 2 

11.0 

22.4 

2.8 

Sodium Iodide 

Nal 

46.6 

111.0 

3.67 

Ilford G-5 Emulsion 

0.4 1H 
0.07226C 
0.01932N 
0.066110 
0.01 189S 
0.349 IBr 
0.474 lAg 
0.003121 

35.8 

81.1 

3.815 


‘Composition listed by percentage weight. 







Table 3. Values of Constants b i and d { 


i 

V 

d* 

i 

1 

0.2335 ± 0.0091 

(2.98 ± 0.30) x 10 3 

2 

1.290 ± 0.015 

6.14 ± 0.29 

3 

(1.78 ± 0.36) x 10 -4 

1.026 ± 0.020 

4 

0.9891 ± 0.0010 

(2.57 ± 0.12) x 10 3 

5 

(3.01 ± 0.35) x 10" 4 

0.34 ± 0.19 

6 

1 .468 ± 0.090 

(1.47 ± 0.19) x 10 3 

7 

( 1.180 ± 0.097) x 10 -2 

0.692 ± 0.039 

8 

1.232 ± 0.067 

0.905 ± 0.031 

9 

0.109 ± 0.017 

0.1874+ 0.0086 


‘From reference 33. 






dn/dco, 1/cm-keV dn/dco, 1/cm-keV 



Electron energy, keV Electron energy, keV 


(a) Spectrum with 0.05 MeV proton in water. (b) Spectrum with 0.1 MeV proton in water. 



Electron energy, keV Electron energy, keV 

(c) Spectrum with 1 MeV proton in water. (d) Spectrum with 10 MeV protons in water. 

Figure 1 . Secondary electron spectrum from incident protons in water. 
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Ion energy, MeV/amu 


Ion energy, MeV/amu 


Figure 4. Electron range as function of ion energy in 
carbon. 


Figure 5. Electron range as function of ion energy in 
silicon. 



Figure 6. Electron range as function of ion energy in gold. 







Figure 8. Theoretical and experimental data of dWtdr in water. 
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(a) For 1 MeV/amu. 



(b) For 10 MeV/amu. 


Figure 9. Effects of electron transmission on calculation of radial dose in water. 
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Dose, Gy Dose, Gy 




(c) For 0.25 MeV/amu 4 He ion. 



Figure 10. Calculations and experimental data of radial dose. 
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Dose, Gy Dose, Gy 



(e) For 2.0 MeV/amu 12 C ions. (f) For 2.57 MeV/amu ,6 0 ions. 




(g) For 90 MeV/amu 56 Fe ions. 


(h) For 377 MeV/amu "°Ne ions. 


Figure 10. Concluded. 
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Dose, Gy 


Radiai distance, ng/cm^ 
(a) For 1 MeV/amu ion. 


Radial distance, ng/cm 
(b) For 10 MeV/amu ions. 



(c) For 1 00 MeV/amu ions. 


Figure 11. Radial dose in carbon, silicon, and gold. 
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